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Abstract 

We study spherically symmetric solutions of the cubic covariant Galileon model in curved space- 
time in presence of a matter source, in the test scalar field approximation. We show that a cos- 
mological time evolution of the Galileon field gives rise to an induced matter-scalar coupling, due 
to the Galileon-graviton kinetic braiding, therefore the solution for the Galileon field is non trivial 
even if the bare matter-scalar coupling constant is set to zero. The local solution crucially depends 
on the asymptotic boundary conditions, and in particular, Minkowski and de Sitter asymptotics 
correspond to different branches of the solution. We study the stability of these solutions, namely, 
the well-posedness of the Cauchy problem and the positivity of energy for scalar and tensor pertur- 
bations, by diagonalizing the kinetic terms of the spin-2 and spin-0 degrees of freedom. In addition, 
we find that in presence of a cosmological time evolution of the Galileon field, its kinetic mixing 
with the graviton leads to a friction force, resulting to efficient damping of scalar perturbations 
within matter. 
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I. INTRODUCTION 



A scalar-tensor theory of gravity with higher-order derivative terms in the action, but 
with no more than second derivatives in the field equations, has first been found in Ref. 
Later, a subclass of the scalar part of this theory has been discussed in a different context 
in j2[. In the framework of the Dvali-Gabadadze-Poratti (DGP) brane model [3|, the cubic 
scalar higher derivative action appears naturally in a specific decoupling limit. More recently, 
a generalization of the DGP decoupling limit action has been presented in |4i], where this 
model was dubbed "Galileon" , because the considered action is invariant under the gahlean 
transformation of the field (adding a constant to the field or a constant vector to its gradient). 
The original Galileon model has been later covariantized to include a dynamical metric |5|, 
and then extended and generalized in several papers, see, e.g. jsl, 0]. 

The fact that the field equations are no more than second order, in spite of the higher- 
order action, guarantees the absence of any Ostrogradsky ghost, i.e., an extra degree of 
freedom with negative kinetic energy. Another interesting feature of the Galileon model is 
that it generically possesses a Vainshtein screening mechanism. This last property can be 
understood as a particular case of "k-moufiage", i.e., of kinetic screening jsj. The Vainshtein 
mechanism hides the propagating scalar degree of freedom, thus providing a way to pass 
local gravity tests, even if the scalar field is directly coupled to matter. The local effects of 
the Vainshtein screening of the Galileon has been studied in a number of works jot] . 

Usually, when studying local effects of the Galileon, the cosmological evolution of the 
scalar field is not taken into account, because of the smallness of the time derivative of 
the field compared to the spacial gradients. However, such a disregard of the Galileon's 



time evolution may be misleading. For example, it was shown in Ref. [10| that in spite 
of the Vainshtein mechanism, which operates locally, this time evolution still poses severe 
constraints on many interesting scalar-tensor theories possessing a shift symmetry and a 
matter-scalar coupling. 

In this work, we study the behavior of the covariant cubic Galileon model (the covariant 
version of the decoupling limit of the DGP model) in presence of a spherically symmetric 
matter source, and possible time evolution of the scalar field imposed by cosmology. We 
investigate in detail two distinct cases: (i) a spherically symmetric matter source is em- 
bedded in asymptotically Minkowski spacetime, and the Galileon field does not depend on 
time; (ii) asymptotically de Sitter Universe with time-dependent Galileon field. To realize 
the latter regime, we adopt the cosmological solution of the Kinetically Gravity Braiding 
model for which the acceleration of the Universe is driven by the time evolving Galileon 
field itself. We find non-homogeneous solutions in the test scalar field approximation, and 
analyze the stability of perturbations on top of these solutions. 

The paper is organized as follows. In Sec. [Tll we introduce the action and field equations 
of the cubic covariant Galileon. In Sec. IIIH we derive a solution for a scalar field having a 
cosmological time evolution, in an arbitrary spherically symmetric static spacetime, in the 
test field approximation. Perturbations are investigated in Sec. IIV[ where we diagonalize 
the kinetic matrix, demixing the spin-0 and spin-2 propagating modes. The solution in 
asymptotically Minkowski spacetime, in presence of a matter source with no time evolution 
of the scalar field, is studied in Sec. IVj where we find the effective metric for the helicity-0 
propagating mode, and the conditions for stability of the solution. In Sec. IVIt we perform 
a similar analysis for asymptotically de Sitter spacetime, in presence of a cosmological evo- 
lution of the scalar field. Scalar perturbations inside a matter source are investigated in 
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Sec. IVIII Besides the analyses of the Cauchy problem and the positivity of energy, we also 
exhibit a large friction term, which is caused by the cosmologically-imposed time evolution 
of the Galileon within matter. Our conclusions are presented in Sec. IVIIII 



II. ACTION AND FIELD EQUATION 

We consider the quadratic plus cubic Galileon action in curved spacetime, 

S = 2Ml j d'x^gi^^- ^{d,vf - ^n^id.^y^ + 5„,atter [V^matter; ~9,.] , (1) 

where the scalar field (p is chosen to be dimensionless, where the reduced Planck mass^ 
Mp = (SttG)"^/^ should not be confused with the mass scale M entering the Galileon's 
cubic kinetic term, and where g and R denote respectively the determinant and the scalar 
curvature of the metric (7^1, (of signature — h ++) used to define covariant derivatives and to 
contract indices. All matter fields, globally denoted as ^/'matter, are assumed to be minimally 
coupled to a physical ("Jordan") metric g^i, = e'^°"^gfj_^, where a denotes a dimensionless 
matter-scalar coupling constant. The constants k2 and k^, dimensionless too, could a priori 
be reabsorbed in the definitions of the scalar field (p and of the mass scale M (while changing 
simultaneously the value of a: y?''"™ = l/tal^/V, M"""^ = M\k2\^/^/\k3\^/^, = a/lfcal^/^). 
In other words, the model actually depends on only two independent parameters, for 
instance M and a. However, the redundant factors k2 and ks will be useful in the following 
to keep track of the origin of the different terms, and above all to change easily the signs of 
these two kinetic terms. Their numerical factors are chosen to simplify the following results 
involving ks, and so that k2 = 1, ks = defines a canonically normalized positive-energy 
spin-0 degree of freedom. Another interest of these coefficients is that any normalization 
convention can easily be recovered, for instance ^2 = 1 corresponding to the relative weight 
of the Einstein-Hilbert action and the scalar-field standard kinetic term in many cosmological 
papers. 



When k^ = 0, i.e., in standard scalar-tensor theories of gravity |12|-|14j|. g^^^ is called the 
"Einstein metric" and its fluctuations describe a spin-2 degree of freedom. As we will see 
in Sec. llVt this is no longer the case when ^3 7^ because of the metric-scalar derivative 
coupling entering the Galileon's cubic kinetic term. The metric (7^1, is thus neither the Jordan 
{g,_iu) nor the Einstein one, and it should just be considered as the variable we choose to 
write the field equations as simply as possible. 

We shall study a test scalar field in the background metric generated by a spherical body, 
assuming that its backreaction on this metric is negligible, which is always the case for a 
small enough matter-scalar coupling constant^. We shall also check a posteriori in which 
conditions this backreaction is indeed negligible. It will thus suffice to focus first on the 
scalar field equation, which can be written as 

V^J^ = -«T, (2) 



^ Throughout this paper, we choose natural units such that h = c = 1. 

^ Anticipating on definition ()10|) below, we are actually talking here about the effective matter-scalar 
coupling constant aeff, which takes into account the induced coupling in presence of a cosmological time 
evolution of the scalar field. 
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where = —(1/ y/^)SS/ Sd^ip is the scalar field's current, and T denotes the trace of T'^'' = 
(2/^/^)(55'niatter/(^5'/ii/), related to the physical (Jordan- frame) matter energy- momentum 
tensor T'^'^ = (2/^^)((55matter/%^) by T^'^ = e^^^T^^ ^ T = e^^^T. The current J'' 
reads explicitly 

^ = 2k,d^^ + - ^ {{d^^f) , (3) 



and Eq. ([2]) takes thus the full form 

k2Uv + ^ - (V,9,(/.)2 - R^'^d.^d^v] = -^T, (4) 

where the Ricci tensor R'^'^ enters because of a difference of third-order covariant derivatives 
of the scalar field: [V^Vj/ — Vj.V^] V^v? = —R'^d^ip. Note that this curvature tensor involves 
second derivatives of the metric, therefore Eq. actually mixes spin-0 and spin-2 degrees of 
freedom. We will diagonalize them when studying perturbations in Sec. HVl but to derive the 
spherically-symmetric solution for a test scalar field, it suffices to replace this Ricci tensor 
by its vacuum value R^^'^ = outside the central massive body. 



III. BACKGROUND SOLUTION 

A static and spherically symmetric background metric reads in Schwarzschild coordinates 

where A(r) and z/(r) are two functions of the radial coordinate. The exterior (vacuum) 
solution is well known to read e'^ = = 1 — r5/r in an asymptotically fiat spacetime, 
where rs = 2Gm denotes the Schwarzschild radius of the central body of mass m, but we 
will also consider the Schwarzschild-de Sitter solution in Sec. |VT] below. We look for a test 
scalar field solution of Eq. 02]) or assuming that the cosmological evolution imposes a 
linear time dependence, i.e., that ip = (where a dot denotes a time derivative). As shown 
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111 ], such a linear time-evolution is a cosmological attractor. Let us look for a solution of 
the form 

ip = (t>{r) + ipct + ipo, (6) 

where ipc and ipo are assumed to be constants. The ansatz has been used in a similar 
context for the study of the Galileon accretion [T^ and the effects of cosmologically evolving 
scalar field on the variation of the Newton constant Actually, since (f' = 0' (where a 



prime denotes a radial derivative), the notation will not be useful in the following, and 
we can thus write Eq. in terms of ip'. With the ansatz this field equation reduces 
in vacuum to a mere dr [r^e*^''*^''^/^ j^-j = g, and it can be integrated once as e^'^^'^^^'^J'^ = 
aMpTs/r'^, the constant of integration being imposed by the same Eq. ([2]) within matter.^ 



Actually, this constant of integration is slightly modified by self-gravity effects within the body, because 
the scalar field is sourced by the trace of the matter energy- momentum tensor, with a different pressure 
dependence than in Einstein's equations. We should thus write rigorously, like in Brans-Dicke theory 
|3,H, e^^+'^^/V = a{l - 2s)Ml,rs/r^, where s = -dX-umjdhiG ~ l-Bgravl/m ~ rs/rtody 
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In terms of this reads 
2h 



v'r 



whose two solutions are 

k2M'^e^r 



ars 



M2 



(7) 



(1 + ^) 



1 ± 
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a 



h .2 

M2 rse 



„ „{X+u)/2 



Note that for a test scalar field, whose backreaction on the metric is negligible, these expres- 
sions dH]) are exact, i.e., correct to any post-Newtonian order (while taking into account the 
slight numerical change of a — a(l — 2s) due to self-gravity, as underlined in footnote |3]) . 
Of course, for the exterior Schwarzschild solution, they can be simplified because e^'^" = 1. 
Moreover, for r 3> rs, one may expand A and u in powers of rg/r, and get 



k2M^r 
4fe 



1 ± 




where 



tteff 



•2 



r 



M2 



(9) 



(10) 



is an effective matter-scalar coupling constant, modifying the bare one, a, because of the 
self-interaction of the scalar field with its own cosmologically-imposed time evolution 0c- 
Let us underline that even if one assumes a = in action ([T]), i.e., a priori no matter-scalar 
interaction, the cosmological time evolution of the scalar field does generate a nonvanishing 
coupling to matter. This is due to the quadratic terms entering the current ([3]), since the 
Christoffel symbols do depend on the background metric and thereby on the mass of the 
central body. 

We will discuss in Sees. |V] and |VT] which sign should be chosen in solution ([8]), by studying 
the large-distance behavior of the scalar field. We will see in particular that it does not reduce 
to ([9]) at cosmologically large distances if spacetime is not assumed to be asymptotically flat 
(see Eqs. (EI]) to fl5^ ). At present, let us just note that the square root entering ^ involves 
a contribution oc 1/r'^, dominating at small enough distances (still assumed to be much 
larger than the Schwarzschild radius rs)- We define the "Vainshtein radius" as 



rv 



tteff 



1/3 



(11) 



It is positive only if k^acs > 0, which is always the case when the bare a = 0, but depends 
on the model when a 0. When ry < 0, then Eqs. (IH]) or (jH]) simply do not give an 
acceptable solution for the scalar field at radii r < |ry|. This means that the ansatz ([6]) 
cannot be correct at such small distances, and the actual solution must mix its time and 
radial dependences in a more subtle way. When ry is positive, on the other hand, Eq. ([8]) 
or ([9]) is a correct solution. In the range ^ r ^ ry, the scalar field then takes the form 



sign{k2k. 




V r / \ry 



(12a) 
(12b) 
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Note in passing that the product /csv?', which will play an important role in the following 
sections, has a sign which does not depend on ^3 (provided this Vainshtein regime exists, 
i.e., notably that k^aes > 0), but which does depend on the sign of /c2, on the contrary. 

It should be underlined that this regime rs r ry can exist even when a = 0, 
provided \2k3(fc/k2M'^\ ^ rs- Therefore, although one could have naively deduced from 
([H])-([nD that if' either vanishes or is proportional to r (depending on the sign to be taken in 
this solution), Eq. fll2l) actually proves it is proportional to r~^/^. The homogeneous scalar 
field assumed in Ref. [11] thus cannot be used in the vicinity of a massive body, when there 
is a nonzero cosmological time evolution 0c 7^ 0. It would remain possible to assume that 
ctes = 0, but this would need to fine-tune the cosmologically-driven time derivative 0c to 
the precise value {—a/k^Y^'^M involving the constants parameters of action ([1]) (while also 
assuming that ak^ < 0). 



IV. PERTURBATIONS 

Let us now study the perturbations of both the metric and the scalar field around the 
background solution defined by Eqs. ([5]) and ([8]). We write g^^^^ = gfj,„ + h^v and ip^^^^ = ip + n, 
where g^u and ip denote the background fields, and h^^ and vr their perturbations. Our aim is 
to identify the pure helicity-2 and degrees of freedom, to analyze in which conditions they 
both carry positive energy, so that any ghost instability is avoided. We expand action ([1]) 
to second order in the perturbations, and keep only the kinetic terms, which involve at least 
two derivatives of these fields. Of course, terms of the form / (background) /iVVvr can be 
integrated by parts, and contribute thus to the kinetic terms as —/ (background) V/iVvr. 
Because of the cubic Galileon interaction in action ([T]), one actually also gets a term 
oc 2dn(fd'^TTnTT, involving three derivatives of the scalar perturbation vr, but this can be 
rewritten as □v9(9^7r)^ — 2V ^dufd^Tid^n by partial integration, as expected because the 
Galileon field equations (perturbed or not) are known to involve at most second derivatives. 
The final sum of all kinetic terms for the perturbations hf^,y and vr reads 



^kinetic 



2M|, 



A;, 



2M2 



2 '7 2 

1 



2 



(13) 



where P"^^^ = U'^^g^^ 



gP'^, and h 



2. 



is the trace of h 



the indices of h. 



fj_y are raised with the background metric 
10 tiic ticn^c yjL it^y. The presence of cross terms / (background) VvrV/i 
illustrates that the helicity-2 and degrees of freedom are mixed. But a simple change of 
variables actually suffices to diagonalize these kinetic terms. We define 



4A;3 



1 



vr, 



(14) 



and Eq. f[T3|l then takes the form 
1 



2Mj, 



/'kinetic 
'-2 



A 

V;A 



(15) 
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where 



->uu 
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2^3 

M2 



2kn k"^ 

^ V^9> + 4^ (d.cpf d^cpd^cp, (16) 



proving that hf^,y describes a pure hehcity-2 field propagating in the curved background of 
Qfj^u, and that the pure hehcity-0 field tt propagates in the effective metric Q'^'^. Our full 
diagonalization recovers thus this effective metric first derived in 11|, from a triangulation 
of the kinetic terms (using the clever trick of replacing the Ricci tensor entering Eq. (jlj) by 
its source, obtained from the Einstein equations). The C(fc|) terms are the only subtleties 
introduced by this diagonalization, but we will see in the following that they are subdominant 
in natural situations. The crucial information brought by Eq. f|T5|) is that the spin-2 graviton 
is never a ghost, and that the effective metric f|T6|) should be of signature — h ++ for the 
scalar perturbation tt to carry positive energy and have a well-posed Cauchy problem. We 
shall perform this analysis in Sees. |V] and \Vl\ for two different asymptotic conditions. 

Before entering this discussion, let us underline that our diagonalization ( |T5i) also allows us 
to exhibit the induced matter-scalar coupling already noticed in Eq. f llUI) for the background 
solution. We can show now that the scalar perturbation vr is also directly coupled to matter, 
even when the bare coupling constant a vanishes. Indeed, since matter is assumed to be 
minimally coupled to ^^"^^ = exp(2a(y9^"")5f|j"" (where "full" means as before the background 
fields plus their perturbations), the action of a point particle reads 



matter 



mc {ffclx^'dx") 



1/2 



, (17) 



where = dx^ / {g^j.dx^dx^Y^'^ denotes the unit 4-velocity of the particle (with respect to 
the background metric (^^j/). Although there does not seem to exist any linear coupling to 
TT when a = 0, we actually know that /i^j, describes a mixing of spin-2 and spin-0 degrees 
of freedom. If we replace it in terms of vr and the actual spin-2 excitation /j,^,^, Eq. (fT^ . we 
thus get 



Sr, 



mc 



a/2 



TT 



(18) 



The quantity within curly brackets plays the role of an effective linear coupling constant 
of matter to the scalar perturbation vr, and it can be nonzero even if a = 0. It should a 
priori not be confused with ttefr, defined in Eq. (fTOj) . which described the effective coupling 
of the background scalar field (f to the matter source. However, at lowest post-Newtonian 
order, they actually coincide when they are constant, notably in a background such that 
0c 7^ but drip = 0. Indeed, we have — {dxipf' = c'^ipl = {u^d\(p) + 0{v), where v is the 
particle's velocity, so that the quantity within curly brackets in reduces to ( TTU]) . 



V. ASYMPTOTIC MINKOWSKI SPACETIME 

We consider in this section that spacetime is asymptotically Minkowskian, i.e., we impose 
that the spherically symmetric metric ([5]) is given by the Schwarzschild solution = = 
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1 — rs/r, and we consistently assume that the background scalar field has no cosmological 
time evolution, ipc = 0. Then the upper sign of solution (E]), ([9]) or (fT2|) must be discarded, 
otherwise the scalar field would diverge at spatial infinity (as well as its derivative and its 
energy-momentum tensor). With the lower sign, Eq. gives for r — )• oo the standard 
behavior of a Brans-Dicke scalar field, 



i + o(— 

\ r 
l + O 



ars 



2k2r 



+ 



V 



(19a) 
(19b) 



In this large-distance regime, the perturbations of the scalar field also behave as in Brans- 
Dicke theory, i.e., they propagate in an effective metric ( |T6l) dominated by its standard k2g'^'^ 
contribution. In conclusion, ^2 rnust be positive for the scalar degree of freedom to carry 
positive energy in the asymptotic r — >■ oo region. 

On the other hand, in the Vainshtein regime r ^ ry (which exists only if k^a > 0), the 
background scalar field takes the form (fT2l) with its lower sign, and the effective metric (fT6l) . 
in which scalar perturbations propagate, is dominated by its k^ contributions. Indeed, we 
have {2k3/M^)Dip ^ |fc2(ry/r)3/2 > > 0. Similarly, the -{2k3/M'^)Vd''f contribution 
to Q^'^ gives at lowest post-Newtonian order ^ k2{rv / r)^^'^ for the Q"' component, minus 
twice this expression (multiplied by g^^ and g'^'^) for the angular Q^^ and Q'^'^ components, 
and again this expression but multiplied by a negligible factor r^/r for the Q^^ component. 
Finally, if we assume that a^//c2 is at most of order one (a^//c2 ^ 1 still being allowed), in 
order to avoid a too large matter-scalar coupling in the asymptotic r — )■ oo region, then the 
assumption rs 'ti r <^ ry implies (/cs^j'^/M^)^ ^ (ar5/4r)^ ^ /c2 <^ {2k3/M'^)D(p, therefore 
the 0{kl) terms entering the effective metric flTB]) are fully negligible: They are of second 
post-Newtonian order beyond the Vainshtein effect. In conclusion, the effective metric (fT6|) 
reads at lowest order 



'Vainshtein 



^2 
T 



r } 



diag —3, 4 



sin^ d 



(20) 



and is has therefore the right — V ++ signature, warranting that the Cauchy problem is well 
posed and that the scalar field carries positive energy. The large numerical factor (ry/r)^/^ 
implies that scalar perturbations are weakly coupled to matter; this is a consequence of the 
Vainshtein mechanism. The different integers entering the diagonal matrix mean that the 
speed of scalar waves is c/ in the orthoradial directions, but 2c/ -\/3 in the radial one. As 
discussed in (igI . 17| . this superluminal radial velocity does not violate causality precisely 
because the Cauchy problem remains well posed simultaneously for all the fields: At any 
spacetime point, there exists a hypersurface which is spacelike with respect to both g^^ and 
Q^^ ^ on which one may specify initial data, and such hypersurfaces may be used to foliate 
the full spacetime. 

To conclude, the curved-spacetime Galileon model ([T]) in an asymptotic Minkowskian 
Universe is consistent both in the Vainshtein regime (r <^ ry) and at large distances, pro- 
vided ^2 > (and fcaa > for the Vainshtein regime to exist). For such an asymptotic 
Minkowski spacetime, the background scalar field is given by solution (IHI), (E]) or (|T2l) with 
their lower sign. 

This conclusion also assumes that o? jk^ is not an extremely large dimensionless number, 
i.e., that matter is not too strongly coupled to the scalar field. If |a| 3> (|/c3|/M^r|)^/^, then 
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the 0{(p'^) contributions dominate in the effective metric (fT6|) at small enough distances, 

and its signature becomes + H instead of — h ++. Scalar perturbations are thus ill 

behaved, not only because their energy can be negative, but above all because their field 
equation is not hyperbolic. However, the present study assumed from the beginning that the 
backreaction of the scalar field on the Schwarzschild metric was negligible, which is obviously 
not the case in the very strong matter-scalar coupling limit. Therefore, this non-hyperbolic 
signature does not even prove the inconsistency of the Galileon model ([T]) in this limit. It 
just means that we must assume as above that a^//c2 is at most of order one, otherwise the 
scalar backreaction on the metric cannot be neglected. 

To estimate this backreaction, it suffices to compare the scalar's energy-momentum tensor 



1 



Ml 



k 



+gfiudp(p V {dxipf - df,ip Vu {d\(pf - dyip {dxff , (21) 

with the one of the matter source, or more specifically their spatial integrals 
{—Tq + T-) 4:7ir''^dr' , which generate the Newtonian potential ^{goo + 1) at a distance 
r from the center of the body. For any r, one then finds that the scalar's contribution is 
O (a^/3(Mrs')^/3) smaller than that of matter. Our analysis can thus be trusted only if a is 
not too large, so that it does not compensate the factor (MrsY^^- On the other hand, even 
if the matter-scalar coupling is of order one (i.e., /k2 ~ 1), then the scalar's backreaction 
is negligible as soon as 1/M is chosen large enough with respect to the Schwarzschild radius 
of any body. In Galileon models, M is generally assumed to be of the order of the Hubble 
constant H (and we will actually derive so in Sec. IVII below, but while assuming a different 
asymptotic Universe). Then Mrs is always an extremely small number, about 10'^^ for 
the Sun, 10~^^ for a galaxy, and still 10~^ for a cluster. Therefore, our test scalar field 
approximation is fully safe if M ~ if and o? Ik2 is at most of order one. 



VI. ASYMPTOTIC DE SITTER SPACETIME 

Let us first consider an isotropic and homogeneous Universe, described by the Friedmann- 
Lemaitre- Robertson- Walker (FLRW) metric ds^ = —dr"^ + a(r)^ {dp^ + p^dVL^), where the 
cosmological time r and the comoving radius p should not be confused with the time t and 
radial coordinate r of the Schwarzschild metric (jS]). Since action ([T]) does not involve any 
cosmological constant, the expansion of the Universe will be caused by the Galileon field ip 
itself. As was shown in [llj, the combined Einstein and scalar field equations indeed admit 
a self-accelerating solution. Let us assume, like in Eq. ([6]) above, that the scalar field has 
a linear time dependence, ip = ipcT + fo, i-e., that (p = (where a dot denotes a derivative 
with respect to r). Then the field equations read 



a{e-3p), (22c) 
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where H = a /a is the Hubble parameter, and e and p denote the energy density and 
pressure of matter, that we shall neglect in the following. [Actually, Eqs. f l22ap and fl22cp 
are valid even if y3 7^ 0.] The only Christoffel contribution to the current (|3]) comes from 
Dip = —SHifc, and Eq. fl22cp reduces to 



which admits the solution (pc = k2M'^ /SHks. [In an expanding Universe, the other solutions, 
= const /a^ — 0, tend either towards this one or towards 0c = 0.] Then Eq. ( ]22bl) 
shows that the scalar-field analogue of e + p vanishes, i.e., that it behaves as an effective 
cosmological constant, and we consistently find that H is also constant. Finally, Eq. fl22ap 
reads 3H^ = — /c20c; implying that k2 must be negative for this self-accelerating solution to 
exist, and it gives the numerical value 



This allows us to rewrite ipc in terms of the constants entering action ([T]). Assuming an 
expanding Universe (i.e., H > 0), we get 



Let us now consider a static and spherically symmetric body embedded in such an ex- 
panding Universe. As before, we assume that the backreaction of the local scalar field on the 
metric is negligible, but its cosmological time evolution ( l25l) is obviously taken into account 
since it is responsible for the accelerated expansion. The local scalar-field solution can still 
be written in the form ([8]), ([9]) or (fT2|) . at least at small enough distances, and the effective 
metric in which scalar perturbations propagate is still given by Eq. (fT6|) . We saw in Sec. |V] 
that the scalar perturbations are ghostlike at large distances if k2 < 0, but this was derived 
while assuming asymptotic flatness and thereby the Brans-Dicke like behavior (fT9|) of the 
background scalar field at infinity. Therefore this previous result is no longer valid in the 
present de Sitter Universe. However, the small-distance physics remains a priori unchanged, 
notably within the Vainshtein radius, and we saw that the signature of the effective metric 
( IT^ depended on the sign of k^dip. As mentioned below Eq. (IT^ . this sign actually does 
not depend on k^, but it does depend on k2, and Eq. (pi]) confirms that we a priori get a 
ghostlike scalar perturbation if ^2 < 0. The stability of the self-accelerating solution seems 
thus spoiled by the presence of any massive body. 

A first way out would be to consider the fine-tuned model such that Ocs = 0. Then 
the scalar field would not be perturbed at all by local matter, there would not exist any 
Vainshtein regime, and the proof in Ref. [ll| that the model is stable would then be valid. 
This assumption Oefj = actually seems more natural in the present self-accelerating Uni- 
verse, since the condition = {—a/k^Y^'^M derived at the end of Sec. IIIII translates as 
a = — sign(/c3)-y/|AJ^]73, involving only constant parameters. However, if we take into ac- 
count the matter sources in Eqs. fl22l) . to derive a more realistic expansion of the Universe, 
then neither H nor ip^ will remain constant, and we will eventually reach an epoch for 
which a is no longer tuned to the right numerical value. Moreover, the self-gravity effects 




(23) 




(24) 




(25) 
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mentioned in footnote [3] mean that the bare couphng constant a needs to be replaced by 
a body-dependent product a(l — 2s), which cannot be fine-tuned to — sign(A;3)y^|fc2j73 for 
all bodies. Therefore, even for this specific model, one expects that the scalar field will be 
directly coupled to matter with a nonvanishing agg, at least at some cosmological epoch, 
and this seems to generically lead to ghost instabilities because 

In fact, the small-distance effective metric fl20l) is not correct in the present expanding 
Universe. Indeed, solution (IT^ does depend on the sign of /c2, but also on the global =F sign. 
We did prove in Sec. |V]that the lower (+) sign gave the correct solution in an asymptotically 
Minkowskian Universe, but this is no longer the case in a de Sitter one. This comes from 
the fact that the FLRW coordinates r and p, useful at cosmologically large distances, do 
not coincide with the t and r coordinates defining the static and spherically-symmetric local 
metric ([5]). In presence of a cosmological constant A = 3iJ^, here mimicked by the self- 
accelerating solution fl2^ . we do know that the exact Schwarzschild-de Sitter solution can 
be written as Eq. ^ with 



= e-^ = \- rs/r - {Hry. 
This can be matched to an asymptotic FLRW coordinate system by defining 



(26) 



P = 



2H 

-Ht 



1 - (Hr) 



2H 



In 



(27) 



Let us also set 



B = 1 



rs 



{He^^pf. 



Then the Schwarzschild-de Sitter metric (IS!)- 0261) becomes 



(28) 



-B 



dT + He'"^pdp 
_ 1 - {He"-pf 
rs 



2Ht 



+ 



+e 



2Ht 



1 + 



rs 



iJ2 [e^rpY 



+ e 
1 

? 
+ 



{dp + Hp dry 



B 



+ p^dn' 



(29a) 



dT^ + 

1 

? 



^rs ^^(l 



H^e^^^p^ 
dp^ + p'^dn^ 



dr dp 



(29b) 



in which one recognizes the asymptotic de Sitter metric ds"^ = —dr"^ + e^^"^ {dp"^ + p^dVL^) 
up to corrections of order 0{rs) caused by the local massive body. Note that we merely 
changed coordinates, and that Eq. ( |29ll still defines the same static solution as (|5l)- (!26l) . in 
spite of the presence of time-dependent exponentials. In particular, the Schwarzschild radius 
still corresponds to the constant e^"^ p = rs, the factor e^'^ compensating the fact that we 
now measure lengths with a varying (comoving) ruler. Note also that the spherical body is 
located at p = at any time, therefore it is comoving (actually static) in the background 
de Sitter Universe. 

To derive the cosmological solution f l2^ -f l25|) . we assumed that the scalar field is homoge- 
neous at large distances in de Sitter coordinates, i.e., that it reads ip = ipcT + (po without any 
comoving radius (p) dependence. This means that in terms of the Schwarzschild coordinates 
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t and r, Eqs. ( 127|1 . ip does depend on r. We explicitly get 



^ct + ^\n [l-{Hrf]+ipo 



Hr 



koM^ 



1 - (Hr) 



3h l-{HrY' 



(30) 



being given by Eq. fl24|) . In other words, the correct solution (|8]) for the background 
scalar field y?', embedded in a self-accelerating Universe, should contain a local r dependence, 
and with the precise factor — A;2M^/3/c3 entering Eq. ( l30l) . Let us prove that this result is 
given by the upper sign of ([H]). Indeed, this solution can no longer be expanded as (jH]) for 
too large radii r, because expressions fl2Bl) for z/(r) and A(r) must now be used, so that 
v'e'^ = —X'e" = Ts/r"^ — 2H^r. We get for such large distances (still assumed smaller than 
l/VSH) 



koM' 



Aks l-UHry 



1 ± 



1 1 - 3{Hr) 
3 1 - (iJr)2 _ 



X 









(31) 



therefore the upper sign recovers exactly the asymptotic behavior derived in Eq. fl30|) from a 
matching with the cosmological solution, whereas the lower sign gives this result divided by 
[2 — 3(ifr)^], i.e., with not only an erroneous factor | by also an incorrect radial dependence. 
Solution ([8]), with its correct upper sign, may now be written for Hr ^ 1 as 




l + O 



(32) 



where the constant — | within the square root comes from a compensation between the rs/r^ 
coefficient entering ([8]) and the large-distance behavior of u'r'^/rs, while using the values ( l24ll 
and ( l25l) for and ipc- Note that this solution is different from Eq. ([9]), which was valid 
for an asymptotically Minkowskian spacetime. In particular, ( l32l) does not reduce to the 
Brans-Dicke solution fll9ap for r ^ ry (while still assuming r ^ ^/H), but gives 



/C2MV 3aesrs\ 

+ — — — I X 



3k. 



2k2r^ 



)x[i + o(!;5) + o(r|)+o(i/V) 



(33) 



with different sign and numerical coefficient for the usual (9(l/r^) contribution, but above 
all involving a much larger 0{r) term, that we already found in Eq. (130!) [it is 0{r^/ry) 
larger than the (9(l/r^) contribution]. 

In the Vainshtein regime (r^ ^ r <^ ry), on the other hand, we still get Eq. ( IT^ . but 
with the crucial information that the upper sign must be chosen: 



" 4A;, 




l + 0(^)+0(J)+0(ifV) 



(34) 



This is the opposite expression as the one we used in Sec. |V] for an asymptotically 
Minkowskian spacetime, and the effective metric in which scalar perturbations propagate, 
at small distances, is thus the opposite of (!20|) [up to negligible corrections of relative order 
0{^rs/r) or 0{\/r^Jr^)^ due to the nonvanishing time derivative (pc we now must take into 
account when computing (|T6l) ]. Since we know that A;2 < for the present self-accelerating 
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Universe to be a solution, this effective metric is therefore again of the correct signature 
— h ++, warranting a well-posed Cauchy problem and the absence of ghosts. 

The signature of this effective metric must also be studied at large distances r ^ ry- 
If one neglects all contributions proportional to the Schwarzschild radius of the massive 
body, Q^^ is most conveniently computed in FLRW coordinates. Then one finds that all 
contributions to Eq. fll6p are of the same order of magnitude 0{k2), including those of the 
usually negligible (dip)'^ terms. The effective metric finally reads Q'^^ = diag(— 2|/c2|, 0, 0, 0) + 
0{rs), meaning that the Cauchy problem remains well posed for scalar perturbations, but 
that their sound velocity vanishes, as was already noticed in [ll]. To avoid any instability, 
the 0{rs) corrections should thus contribute positively to the effective spatial metric Q^^ . 
We thus need to compute them, while taking into account the zeroth-order terms 0{r^) 
although we know they will eventually cancel. Indeed some of these zeroth-order terms get 
multiplied by when solution ([H]) is expanded to compute flTBl) . Another complication is 
that Q^^' 7^ in Schwarzschild coordinates fl^-f l2B]) . i.e., that Q'^'^ is not diagonal. However, 
to determine its signature, it suffices to consider the quadratic form Q^ydx^dx^ (where the 
effective metric's indices are lowered with g^y, i.e., Q^y is not the inverse of Q^^), and 
to note that Q^odt^ + 2QQrdtdr + Qrrdr'^ = Goo (dt + Gordr/GoQ) + {Grr — Gor/Goo) dr"^- In 
this new frame diagonalizing the effective metric, one finally finds at lowest order G^iu ~ 

diag (— 2|A;2|, 2|A;2|e, — |/c2|e'"^5 — |^2|e'"^ sin^ 6') , where e = |(rv/r)^, i.e., a signature — I 

instead of — h -|-+, implying that the Cauchy problem is ill-posed for scalar perturbations 
propagating in the orthoradial directions (or that they are unstable if one interprets these 
negative signs as an imaginary sound speed, such instabilities being large at moderately large 
distances of a few ry). In conclusion, although we found above that the physics of scalar 
perturbations is consistent in the Vainshtein regime r ^ ry, it is not at large distances 
r 3> ry- 

This serious problem can a priori have several solutions. The first naive idea would 
be to choose the theory parameters such that ry > 1/-^, so that any instability would 

1/3 

be hidden behind the cosmological horizon. But since ry oc , this cannot work for all 
massive bodies, and this would need anyway a very large matter-scalar coupling constant a, 
thereby spoiling the hyperbolicity in the Vainshtein regime for the same reasons as at the 
end of Sec. |V] (and being inconsistent with our assumption of a negligible backreaction of 
the local scalar field on the background metric). A much better argument is that the scalar 
field may never get out of the Vainshtein regime if one takes into account the full matter 
distribution in the Universe: At any location, there would be a close enough and massive 
enough body such that r < ry. But the strongest argument is that the vanishing sound 
speed found above at order C(r°) is a consequence of the exact self-accelerating solution 
As soon as one takes into account matter in the field equations (12^ . then the 
sound speed becomes nonvanishing and real, as was shown in and it dominates at large 
distances over any 0{rs) correction. For instance, if we assume a = like in the 



presence of a positive matter density e in Eq. ( I22ap increases the value of H with respect 
to the lowest-order expression f l24|) . This implies a modification of Eq. f l25|) for 0c, and the 
matter field equation e + ?>He = (assuming a vanishing pressure p) allows us to compute 
(fc- We finally find that the terms proportional to dip and {d\(p)'^ in Eq. f|T6|) both give 
positive 0{e) contributions to the spatial effective metric G^-' , and we explicitly obtain 





€ 



(35) 
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Surprizingly, the result is fully different if we assume a nonzero bare coupling constant a, 
because the matter field equation e + "iHe = aeipc implies a different time evolution for e, 
and the time integration of Eq. f l22cp makes a go away. One finds that the lowest-order 
expressions of H and (pc both get multiplied by the same factor 1 — e/12MpH^ + 0{e'^), and 
the matter field equation can again be used to compute ipc- Now the term proportional to 
Dip in Eq. (IT^ happens to give a negative 0{e) contribution to Q''^ , but it is counterbalanced 
by the {dxip)"^ and W^d^Lp terms of (ITBi) . We finally get 



which is positive definite if 3a^/|A;2| < 1 (and for even larger |a| if k^a > 0). The reason why 
does not tend to when a — )• is because they correspond to different cosmological 
initial conditions. In any case, their positivity shows that the ill-posed Cauchy problem 
found above was the consequence of an oversimplified cosmological background, and the 
signature of the effective metric at large distances is in fact — h ++, as needed. 

To conclude, the curved-spacetime Galileon model ([1]) with k2 < admits a self- 
accelerating cosmological solution, and a spherical body embedded in this Universe generates 
a scalar field given by the upper sign of Eq. (jS]), i.e., by fl52]) . fl55]) or f l51|) depending on the 
distance to this body — the Vainshtein regime 0341) existing only if k^aeg > 0. Perturba- 
tions around this solution carry positive energy and have a well-posed Cauchy problem in 
the Vainshtein regime r ^ ry. Farther away, one must take into account the matter content 
of the Universe to derive the cosmological expansion, and the effective metric in which scalar 
perturbations propagate then remains of the correct hyperbolic signature — h ++. 

Let us finally mention in which conditions the local scalar field's backreaction on the 
metric is indeed negligible, as was assumed to draw the above conclusions. First, this 
is always the case if cteff is small enough, i.e., if one tunes the bare matter-scalar coupling 
constant to a ~ — sign(/i;3) a/| /c2 1 /3, so that it almost compensates the cosmologically induced 
coupling k^ipl/M'^ entering Eq. (ITU]) . Indeed, the spatial integral of Eq. (1^ over a sphere 
of radius r < is at most of order 0{aesrs), negligible with respect to the effect of 
the material body of mass ^r^. However, for a random bare coupling constant such that 
Q^^/|^2| ~ 1 (a small value being also allowed), Eqs. ffTOj) and fl25|) imply that alf^/\k2\ is 
generically of the order of unity, therefore one actually expects large deviations from our 
results above. But even in such a case, the scalar field's backreaction is still negligible in the 
Vainshtein regime r <C ry, because it is of order O (aesrsir/rv)^^'^) ■ Therefore, our previous 
conclusion that the Galileon model ([T]) is stable in the Vainshtein regime remains valid even 
for a matter-scalar coupling of order one. On the other hand, for such a 0{1) coupling, our 
test scalar field approximation cannot be trusted at distances r ^ ry, and further work is 
needed to actually prove stability jisj. Finally, for cosmologically large distances r ~ 3> 
ry, both the effects of the massive body and the locally generated scalar field are negligible, 
and we asymptotically reach an effective metric Q^'^ ~ \k2\dia:g{— 2, 0{1) e g^^ /MpH"^), with 
the correct hyperbolic signature — h ++. 

VII. SCALAR WAVES WITHIN MATTER 

Although we have shown in Sees. IVland lVII that the field equations for scalar perturbations 
are hyperbolic in the Vainshtein regime, and thereby that there is no ghost instability at 




(36) 
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small enough distances from a massive body, the growth of such perturbations within the 
body deserves a more careful study. First of all, the vacuum solution ([8]) for the scalar field 
is no longer valid inside matter, therefore the hyperbolicity of the effective metric (fT6l) needs 
to be checked. To sirnplify, let us consider a spherical body of constant density. As derived 
in page 331 of Ref. [l9|, the interior solution for the metric (neglecting cosmological 
corrections of order H'^r'^) reads 



.A(r) 



1 

4 



3(1-^ 

2\ -1 



1/2 



- 1 



2\ 1/2' 



n 2 



l + ^ + 0(r|) 



(37b) 



where r* denotes the body's radius, not to be confused with its Schwarzschild radius rs- The 
background scalar field equation ([2]) can then be solved as in Sec. IIIIl with the difference 
that only the mass interior to the sphere of radius r, m(r) = {rs/2G){r/r^)^, is a source for 
(p'{r). We get 



interior 



4fe 



1± Wi- 



^1 + 



aA+O (r|) + O {H'r') , (38) 



where either H = in asymptotic Minkowski spacetime, or H and ipc are given by Eqs. (124 
and (125|) in asymptotic de Sitter spacetime (yielding the constant — | we already encountered 
in Eq. ( 132|) for the exterior solution). The only difference with respect to the exterior solution 
([9]) or (!32|) . at this lowest post-Newtonian order, is that ctefr is multiplied by the constant 
1/rl instead of 1/r^. In the Vainshtein regime r < <C ry, we have thus 



interior 



(?) 



3/2 



(39) 



instead of Eq. (fT2|) . where the lower sign (+) corresponds to asymptotic Minkowski space- 
time, and the upper one (— ) to asymptotic de Sitter spacetime. The various contributions 
to the effective metric f|T6l) can now be computed within the body, and one finds at lowest 
order (in both asymptotic Minkowski and de Sitter cases) 



'interior 



\k. 



(^) 



3/2 



diag 



1 



2 ' ' ' sin^ 9 



(40) 



instead of ( l20l) [or the opposite of of ( l20i) in the asymptotic de Sitter case, for which —k2 = 
\k2\]- Note that this effective metric is discontinuous at the surface of the body, because it 
involves ip", which is discontinuous when one suddenly passes from empty space to a finite 
matter density. But of course, everything becomes continuous when one considers a smooth 
matter profile. Note also that the sound speed, a/2/3 c, is now isotropic, and differs from its 
value outside matter (that we found in Sec. |V]to be 2c/ -\/3 in the radial direction and c/ a/3 
in the orthoradial ones). This can be interpreted as different refractive indices of matter 
and vacuum for scalar waves. But the crucial information brought by Eq. f HOj) is that the 
effective metric remains of hyperbolic signature — h ++ within matter, therefore no ghost 
instability develops even in the interior of a body. 
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On the other hand, the Ricci tensor entering Eq. (jl]) no longer vanishes within matter, 
and this causes extra couphngs of the perturbations to the background fields, whose effects 
could a priori lead to other types of instabilities. This Ricci tensor is actually generated 
both by matter and by the background scalar field itself. The scalar-scalar self-interactions 



are responsible for the 0{kl) terms in the effective metric f ll6p . as was first shown in 11 
Not only we already took them into account to derive fl40l) . but we also saw in Sees. IVl and 
|Vl]that they are negligible with respect to the 0(^3) terms, provided a'^/\k2\ is at most of 
order 1 (i.e., that we are not considering a extremely large matter-scalar coupling constant). 
We can thus focus on the Ricci tensor generated by matter, -Rmatter = | ('"5/'"*) 9^'^ + ('^Di 
which multiplies first derivatives of the scalar field in Eq. and thereby does not change 
its kinetic term (second derivatives). Nevertheless, its presence means that within matter, 
scalar perturbations acquire a direct derivative coupling to the background scalar field. The 
first-order expansion of (jl]) reads 

Dn - (V^9V) (V^S.vr) - RZttcr ¥^,,^,u = 0, (41) 

where we have neglected the 0{kl) terms coming from the diagonalization (fT4l) -ffT6l) of the 
kinetic terms, as well as the {k2M'^ /2k^) nvr linear term because we assume the star is in 
the Vainshtein regime (r* <^ ry). At lowest post-Newtonian order, the radial contribution 
"-^matter V''^' "^^^ ^^^o bc ueglcctcd, because both -R^atter — ^ ('^s) ip' = O (\/rs) tend 
to with rs- But the temporal contribution — -Rmatter 0c ^ is crucial, as it corresponds to a 
damping or antidamping term depending on its sign. Let us thus consider the asymptotic 
de Sitter case of Sec. |VT1 where ipc is given by Eq. (125|) . The main contributions to Eq. ( BTj) 
then read 

(&rior^^^) - ^matter 0c 71" = 

O / I ; I \ 1/4 / \ 1/2 



The coefficient of the tt term is generically of order \/rsJrf, notably when the bare matter- 
scalar coupling constant a vanishes or is negligibly small. Indeed, Eqs. ( ITO|) and (125|1 then 
give ttgfj ~ I ^2 1/3. Note that this coefficient can never be small with respect to \/rs/rl, 
because we assume that a^/|/c2| is at most of order 1, therefore agfr/l^sl is never large either. 
On the other hand, it is possible to choose the bare a so that is small (see Sec. IVip . and 
the 77 term of Eq. ( H2i) may thus be multiplied by a large number in some specific situations. 
Therefore, the extra coupling of the perturbation vr to the cosmologically-imposed ipc within 
matter can change drastically its behavior. The good news is that its sign implies this is a 
friction term, and therefore that no scalar instability occurs within the material body. The 
plane- wave solutions of Eq. P2l) read indeed (in Cartesian coordinates) 



TT oc e sin {ut — k ■ x + const) , (43) 



with a dispersion relation 
and a decay time 



^^ = -^ + ^> (44) 
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For the Sun, 2\Jr^/rs is about one hour, therefore this damping is always quite efficient. But 
when tteff happens to be small, because of a balance between the bare matter-scalar coupling 
constant a the cosmologically-induced one k^Lpl/M"^ in Eq. ffTOj) . scalar perturbations are 
suppressed even faster. For a; = and | < there exist other solutions of the 

form TT oc e~*/^* sin (k ■ x + const), where 1/T± = 1/T ± a/I/T^ — 2^? /?> is always positive, 
therefore such perturbations are damped too. The upper sign gives an even faster decay 
than for the plane waves ( l43l) . whereas the lower sign corresponds to a slower damping. The 
limiting case of a homogeneous perturbation (k = 0) gives either a fast damping vr oc e~^*/^, 
or an actual constant which can be reabsorbed in the definition of the background (p^. 

In conclusion, the physics of scalar perturbations is fully safe in the interior of a spherical 
body in the Vainshtein regime. Not only the Cauchy problem is well posed and the kinetic 
energy carried by scalar perturbations is positive, but even the subtle derivative coupling 
caused by the nonvanishing Ricci tensor in Eq. (HTl) efficiently damps these perturbations. 
Note that this local friction within material bodies is paradoxically caused by the slow 
cosmological time evolution (0^ ^ 0) of the background scalar field. The large-distance 
behavior of the Universe is thus significantly influencing local physics. 

VIII. CONCLUSIONS 

In this paper we studied in detail spherically symmetric solutions for the cubic covariant 
Galileon model ([T]) in presence of a matter source. An important ingredient of our study is 
that we take into account the time variation of the scalar field, induced by the cosmological 
evolution. Although cosmological time evolution is tiny (normally of order of the Hubble 
rate, (pc ~ H), nevertheless its effect on the behavior of the scalar is crucial and cannot be 
disregarded. 

We exhibited, in particular, an «n(iMce(i matter-scalar coupling for (pc 7^ 0, Eq. f lTO|) . This 
effect arises because of the kinetic mixing of the Galileon and the graviton, which results 
in the coupling of the derivatives of the scalar field to the curvature, in the Galileon field 
equation. For (p^ ^ 0, this term plays the role of a matter-scalar coupling, hence it is induced 
by the cosmological evolution of ip. It is important to note that this coupling is naturally of 
order of unity, therefore the Galileon effectively couples to matter sources even if the bare 
scalar- matter coupling is absent, a = 0. 

We also found that the local solution for the Galileon field crucially depends on the 
asymptotic conditions at large distances. Indeed, to find time-dependent solutions in pres- 
ence of massive sources, we applied the ansatz ([6]), which enables to separate the time and 
radial variables. The radial- dependent part can then be found exactly, in the test scalar field 
approximation, for an arbitrary static background metric, Eq. ([H]). [We checked a posteriori 
the conditions for backreaction to be negligible.] The full solution for however, consists 
of two different branches. The solution with the lower sign in ([H]) corresponds to the asymp- 
totically Minkowski spacetime, with no time evolution of ip. The other case we studied in 
detail is the asymptotically de Sitter Universe, whose self-acceleration is provided by a time 
evolution of the Galileon itself. In this case, the cosmological asymptotic behavior dictates 
to choose the other branch of the solution, namely, Eq. ([8]) with the upper sign. It turns out 
that the stability of the solutions drastically depends of the choice of the branch. 

In order to investigate stability of the found solutions, notably the well-posedness of 
the Cauchy problem and the positivity of energy, we identified the actual spin-0 and spin- 
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2 degrees of freedom by diagonalizing their kinetic terms in Sec. IIVI Using the results of 
Sec. HVt we found in Sec. |V]that the background scalar field given by solution ([8]), ([9]) or 
(fT2|) with their lower sign (corresponding to the asymptotically Minkowski spacetime), is 
consistent both in the Vainshtein regime (r ^ ry) and at large distances, provided ^2 > 
and k^a > 0. 

The asymptotically de Sitter case has been treated separately in Sec. |Vll It should be 
noted that the self-accelerating de Sitter stage is only possible for the choice /c2 < in action 
([1]). Therefore naively one would expect appearance of a ghost instability at least in the 
Vainshtein regime, see e.g. Eq. fl20|) . where a negative /c2 means ghost scalar perturbation. 
One should also have in mind that the existence of the Vainshtein regime is generic, even 
with no direct scalar-matter coupling, because of the induced coupling for time-dependent 
solutions. It turns out, however, that the choice of the correct branch, corresponding to the 
asymptotically de Sitter solution, namely, Eq. (jS]) with the upper sign (in contrast to the 
asymptotically Minkowski where the lower sign is taken), changes drastically the behavior 
of perturbations. Indeed, we found that the solution is stable for the Vainshtein regime, 
contrary to naive expectations. 

The large distance behavior, r 3> ry, for the asymptotic de Sitter is more subtle. The pure 
cosmological perturbations of the spin-0 degree of freedom have positive-energy dust-like 
behavior (i.e., with vanishing sound speed). Therefore any small deviation can potentially 
create a change in signature for perturbations. In fact, the deviations caused by a local 
material body do not yield a correct signature: The field equations for perturbations become 
elliptic in the angular directions. However, this problem exists only in the pure de Sitter 
universe filled with the Galileon field with no external matter. In a more realistic scenario, 
when a small amount of normal matter is added (such that the universe is almost de Sitter), 
the effective metric f|T6l) always has the correct signature. 

Another interesting effect, associated with the cosmological time evolution of the Galileon, 
is the appearance of a significant friction term in the equation of motion. Because of the 
commutation of covariant derivatives, the Galileon field equations involve couplings of the 
scalar field to the curvature tensor. In the cubic model ([1]) studied in this paper, this gener- 
ates within matter an extra interaction of scalar perturbations with their background. We 
showed that it causes an efficient damping of these perturbations, a local effect paradoxically 
caused by the cosmological time evolution. 

To summarize, we studied spherically symmetric solutions of the cubic covariant Galileon 
model in presence of a matter source and a cosmological evolution. We found that the 
time evolution of the Galileon, even as small as a cosmological one, leads to considerable 
effects: in particular, the appearance of an induced coupling, which is generically of the 
order of unity, so that the Galileon becomes effectively coupled to local matter sources even 
if the bare coupling is zero; and the emergence of a friction term, which effectively damps 
perturbations within matter sources. The detailed analysis of perturbations showed that 
the Galileon model is well behaved in asymptotically Minkowski space as well as in the 
asymptotically de Sitter, provided that the correct branch of the solution is chosen. 
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